Effects of spin non-collinearities in magnetic nanoparticles 



H. KachkachE 

Groupe d'Etude de la Matiere C'ondensee, Universite de Versailles St. 

45 av. des Etats-Unis, 78035 Versailles, France 



Quentin, 



In a many-spin approach that takes account of the internal structure, microscopic interactions and 
single-site anisotropics, we investigate the effect of spin non-collinearities induced by the boundary 
and surface anisotropy on the behaviour of individual magnetic nanoparticles. Through analytical 
and numerical calculations, we show that there are mainly two regimes separated by some critical 
value of the surface anisotropy constant Ks which controls the intensity of spin non-collinearities: i) 
the so called Stoner-Wohlfarth or Neel-Brown regime of a macrospin undergoing a coherent switch- 
ing, ii) the many-spin regime where the strong spin non-collinearities invalidate the coherent mech- 
anism, and where the particle's magnetic state and switching mechanisms can no longer be modeled 
by a macrospin. For small-to- intermediate values of Ks, and within two models of surface anisotropy 
(transverse and Neel), the behaviour of the nanoparticle can be modeled by that of a macrospin 
with an effective potential energy containing a uniaxial and cubic anisotropy terms. This effective 
spin model provides a crossover between the two regimes. 
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I. INTRODUCTION 

Reducing the size of magnetic nanoparticles in view of 
room-temperature technological applications costs a high 
price. Very small particles present two major difhculties: 

i) they are less stable against thermal fluctuations, and 

ii) they have an important contribution from the bound- 
ary which induces large deviations from the homogeneous 
(coUinear) magnetic state. As an immediate consequence 
of the latter the magnetic state of the particle cannot be 
represented by a macroscopic magnetic moment (single- 
moment approach), and one has to resort to finer (local) 
approaches involving the atomic magnetic moment right 
from the beginning. Only then can one distinguish, and 
thereby estimate the role of, the various environments in- 
side of the particle, and in particular one may assess the 
effect of the boundary/surface contribution on the static 
as well as dynamic behaviour of the magnetic particle. 
This of course assumes that one has adequately adapted 
the computing methods. 

The magnetic structure and switching mechanisms of 
many-spin particles (MSP) have been studied for a few 
years now by many authors using numerical methods, 
such as the solution of Landau-Lifshitz equation and 
Monte Carlo technique 0| ■ Analytical calculations using 
the modified spin-wave theory 01 and the spherical model 
Q have shown that the spin disorder on the boundary 
is of long range and propagates deep into the particle. 
Most of the results emphasize the important role of sur- 
face anisotropy and boundary effects and the necessity 
to take account of the internal structure of nanoparticles. 
However, the study of the dynamics of nanoparticles in 
the many-spin approach presents tremendous difficulties 
related with the analysis of the energyscape (minima, 
maxima, and saddle points), which is a crucial step in 
the calculation of the relaxation rate and investigation 



of the magnetization reversal at finite temperature. One 
may then address the question as to whether there ex- 
ist some cases in which the full-fledged theory that has 
been developed for the one-spin-particle (OSP) approach 
[sec and references therein! can still be used to de- 
scribe an MSP. In Refs. 0, 0, II we showed analytically 
and numerically that indeed, for small deviations from 
the collinear state, we can model the magnetic state and 
switching mechanisms of an MSP by a macrospin in an 
effective potential energy containing a uniaxial and cubic 
anisotropy terms. 

In this paper, we will briefly review through some ex- 
amples that results that show how the surface defects, 
and the spin non-collinearities they entail, fundamen- 
tally alter the magnetic properties of a small (3 — 10 nm) 
nanoparticle. In particular, we will show that such a 
nanoparticle behaves according to different regimes with 
crossovers determined by the underlying material param- 
eters, notably the surface anisotropy (model and inten- 
sity). It turns out that in some typical parameters ranges, 
namely weak surface anisotropy, the OSP approach may 
be rehabilitated provided that an effective picture of the 
particle is adopted. More precisely, the MSP may be 
represented by a single (macroscopic) magnetic moment 
with an effective energy. This simplification provides us 
with a useful tool for studying the dynamics, and the 
thermally activated switching of the magnetization, of 
a small nanoparticle while taking account, though in a 
phenomenological manner, the effect of its boundary. 



II. MODEL AND NOTATION 

Surface effects are due to the breaking of crystal-field 
symmetry at the boundary of the nanoparticle. In or- 
der to study such effects, one has to resort to micro- 
scopic theories capable of distinguishing between differ- 
ent atomic environments and taking account of physi- 



cal parameters, such as surface anisotropy, exchange and 
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magneto-crystalline anisotropy in the core and magnetic 
field. 

Our targeted system is an MSP of M (reduced) atomic 
magnetic moments mi with |mi| = 1. The particle's en- 
ergy includes the exchange, the Zeeman, and anisotropy 
contributions 



H 



■E 



'H.-mi + KiA{m.,) + Y^ 



A: 



a) 

where A(si) is the anisotropy function whose expression 
depends on the type of anisotropy, it may be, e.g., of the 
uniaxial, cubic, or Neel's type, that is 
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with u.y being a unit vector connecting the nearest neigh- 
bors z, J. The anisotropy constant Ki may be positive or 
negative and is denoted by if the site i is in the core 
of the particle and if it is on the boundary. A spin in 
the core has its full coordination number while a spin on 
the boundary lacks some of its neighbors. 

In the present work, the anisotropy in the core is taken 
as uniaxial with the anisotropy axis along the reference 
z axis. For surface spins, we take either the (uniaxial) 
Transverse surface anisotropy (TSA) model with an axis 
along the radial (i.e., transverse to the cluster surface) 
direction or the Neel surface anisotropy (NSA) model 

1^ ^3 ■ We use the more general model of transverse 
direction given by the gradient [the vector perpendicular 
to the isotimic surface ?/; = constant defining the shape 
of the particle, e.g. a sphere or an ellipsoid]. In the case 
of a spherical particle, the transverse and radial direc- 
tions coincide, whereas for another geometry such as an 
ellipsoid they do not [see Ref. for more details]. 

In the sequel, we will make use of the dimensionless 
parameters = KclJ.,ks = Ks/J,t = kBT/J,h = 



III. MANY-SPIN PARTICLES 

A. Effects of non-collinearities 

The l/D surface contribution to AV,off, where D is 
the particle diameter, is in accord with the picture of 
all magnetic atoms tightly bound by the exchange in- 
teraction, whereas only the surface atoms experience the 
surface anisotropy. This is definitely true for magnetic 
films where a huge surface contribution to the effec- 
tive anisotropy has been observed. The same holds for 
cobalt nanoclusters of the form of truncated octahedrons 
|lOj . where contributions from different faces, edges, and 
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FIG. 1: Magnetic structure of a spherical nanoparticle of 
linear size A'' = 15 with ks = 2 for the global magnetization 
directed along [1,1,0], showing atoms in the plane z = 0. 



apices compete resulting in a nonzero, although signifi- 
cantly reduced, surface contribution to Ky^cS- However, 
for symmetric particle shapes such as cubes or spheres, 
the symmetry cancels this (first-order) contribution. In 
this case, one has to take into account deviations from 
the collinear spin state that result from the competition 
between the surface anisotropy and the exchange inter- 
action J [see Fig. ^for the NSA]. In the case Kg ^ J 
deviations from coUinearity are very strong, and it is dif- 
ficult, if not impossible, to characterize the particle by a 
global magnetization suitable for the definition of the ef- 
fective anisotropy. On the other hand, in the typical case 
Ks ^ J the magnetic structure is nearly collinear with 
small deviations that can be computed perturbatively in 
Ks/J ^ 1. The global magnetization vector mg can be 
used to define the anisotropic energy of the whole par- 
ticle. The key point is that deviations from coUinearity 
and thus the energies of the system are different for dif- 
ferent orientations of mo, even for a particle of a spherical 
shape, due to the crystal lattice. 
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FIG. 2: Critical field against the ratio of surface anisotropy 
constant to exchange, ka for a spherical particle of 360 spins. 
The anisotropy is uniaxial in the core (with fee = 0.01) and 
TSA on the surface. 

In Fig. 121 wc plot the critical field as a function of the 
ratio of the reduced surface anisotropy to the exchange 
coupling for a spherical particle of 360 spins with uniax- 
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FIG. 3: Magnetization as a function of field with variable 
temperature (t) for a particle of 515 spins with core uniaxial 
anisotropy (kc = 0.0024) and transverse surface anisotropy 
{ks = 0.4). 



compute the contribution of surface anisotropy in Eq. 
we replace the number of nearest neighbors of a surface 
atom by its average value 

Zia => ZiQ = 2 - \na\/ max{|n^|, \ny\, |n^|} . (3) 

where ria is the a-component of the normal to the surface 
n. The surface-energy density can then be obtained by 
multiplying the contribution in Eq. by the surface 
atomic density /(n) = maxjln^;!, \ny\, \nz\}: 

Es{m,n) = -Ks{m-nf fin) (4) 

At equilibrium, in the continuous approximation the 
Landau-Lifshitz equation reduces to 

m X Hoff = 0, Heff = B.A + JAm, (5) 

where A is the Laplace operator and the anisotropy field 
Ha contains contributions from the core and the surface 



ial anisotropy in the core and TSA on the surface. Each 
point on this curve is obtained from the hysteresis loop 
computed for the corresponding value of kg (see Ref. [l| 
for the details of the procedure). These results show 
that upon varying the surface anisotropy one observes 
a crossover from i) the coherent-reversal regime assum- 
ing the particle as a macrospin according to the model of 
Stoncr-Wohlfarth, into ii) the incoherent-reversal regime 
with clustcr-wisc switching. In the latter regime, the 
particle exhibits (due to strong surface anisotropy) new 
features that are reminiscent of a many-spin system, 
which cannot be described by a macroscopic approach. 
In Refs. Q it was shown that for the TSA and NSA there 
exists a (different) critical value of the surface anisotropy 
constant that separates the above mentioned regimes. 

Fig-Elsliows the magnetization of a many-spin particle 
as a function of magnetic field at different temperatures. 
These results, obtained by Monte Carlo simulations (see 
Ref. for details), simply show that, because of the sur- 
face effects, the magnetization does not saturate even at 
relatively high fields, i.e. h = 0.2 which corresponds to 

~ 16 Tesla for cobalt particles even at very low tem- 
perature, i.e. t = 0.01 (T ~ 1 K). 



B. Crossover to an effective one-spin particle 

With the desire to avoid the difficulties mentioned 
above, one may address the question as to what ex- 
tent the dynamics theory developed in the OSP ap- 
proach can still be applied to an MSP. However, avoid- 
ing somehow the spin non-coUincaritics induced by sur- 
face/interface anisotropy means that some price has to 
be paid. 

Let us now briefly summarize the results obtained in 
this respect. Consider a spherical particle with uniaxial 
anisotropy in the core and TSA on the surface. Then, to 
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where N is the side of the cube inside which the particle 
is cut. 

For Kg <C J the deviations of m(r) from the homo- 
geneous state mo are small and one can linearize the 
problem as follows 

m(r) = mo + -0(r, mo) mo + ei-0i 4- e2ip2, 

V = 1^1 « 1. (7) 

The correction i/' is the solution of the Hclmholtz equa- 
tion with boundary conditions 
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where n = r/R. 

The solution -0 of Eq. (jSJ has the form |^ 

V'(r,m) = — / <fr'G{Yy)i{xn,n) 
47r Js 



(9) 



where G'(r, r') is the Green function of the problem. 

In the absence of core anisotropy, where the Hclmholtz 
equation reduces to a Laplace equation, an extact ex- 
pression (call it G(°)(r,r')) for G(r,r') was obtained in 
Ref. . In the presence of core anisotropy we have 
recently found an approximate expression for the 
Green's function which contains a correction to G(r, r') 
of the order oc K^/J ^ 1. This correction is then 
rewritten as a convolution of two G^''-''s. Collecting all 
contributions from the core and surface, we showed that 
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the effective energy of an MSP particle, in the absence of 
magnetic field, is written as 

£efi. = ml + Ki{ml + ml + m^). (10) 

The coefficient K2 of the second-order contribution is in 
fact the result of two contributions, one coming from the 
initial core uniaxial anisotropy and a new contribution 
that is induced by the surface anisotropy. The latter 
contribution is much smaller than the former because its 
coefficient contains the product [Kc/ J){K1 / J) <^ 1. The 
details of the calculations will be published elsewhere Q . 
The 4*''-order coefficient ^4 was found in Ref. to be 
given by 



where Af, Kg, z, J are respectively the number of atoms, 
the surface anisotropy constant (transverse or Neel) , the 
coordination number, and the exchange coupling of the 
many-spin particle, k is a surface integral that depends 
on the underlying lattice, the shape, and the size of the 
particle and also on the surface-anisotropy model. For 
a spherical particle (of ~ 1500 spins) cut from a simple 
cubic lattice with Neel's surface anisotropy, k ~ 0.53465. 

In we confirmed this result by numerical calcula- 
tions of the field behaviour of the net magnetization and 
effective energyscape. Before we discuss the results, we 
briefly explain the method we used to obtain them. Be- 
cause we are dealing with an MSP, the energyscape can- 
not be represented in terms of the coordinates of all spins. 
Instead, we may represent it in terms of the coordinates 
of the particle's net magnetization. For this purpose, we 
fix the global or net magnetization, m, of the particle in 
a desired direction mg (|ttio| = 1) by using the energy 
function with a Lagrange multiplier A 0: 

V s 

T = n-AfX-{m-mo), m = (12) 

To minimize we solve the evolution equations 

Si = - [sj X [si X Fi]] , Fi = -dT/dsi 

A = dJ'/dX^-Nim-mo), (13) 

starting from = rriQ = m and A = 0, until a stationary 
state is reached. In this state m = niQ and [s^ x F^] = 0, 
i.e., the torque due to the term 7VA- (m — mo) in T com- 
pensates for the torque acting to rotate the global mag- 
netization towards the minimum-energy directions [see 
discussion in Ref. Q]. The orientation of the net mag- 
netization is then given either in Cartesian coordinates 
(m^:, TOj^, 771^) or in spherical coordinates (0„, </?„). 

We now discuss the main results. The first panel of 
Fig. El (fcs = 0.1) shows that, for very small kg, the en- 
ergyscape of an MSP is well recovered by the effective 
energy in Eq. (|10|) . As kg increases [see middle panel, 
kg = 0.3], some deviations start to be seen, and for rela- 
tively larger values of kg a fit with Eq. (|10|l is no longer 



possible. In fact, in this regime strong deviations from 
coUinearity develop, especially near maxima and saddle 
points, as can be seen on the panel with kg = 0.5 in Fig.E] 
In fact, in this case the Lagrange-parameter method in- 
troduced in 1^ fails because the magnetic state of an 
MSP can no longer be represented by a net magneti- 
zation. These results imply that the effect of the spin 
non-coUinearities on the energy is to split the minimum 
at 6'n = 0, defined by the uniaxial anisotropy, into four 
minima at 9n ~ 28° and ipn — 0, ±7r/2, tt, reminiscent 
of cubic anisotropy [see Fig. |3] . These minima are con- 
nected by saddle points at ipn = ±7r/4 and ±37r/4 and 
the point at 0„ = becomes a small local maximum. 
The four minima exist over a finite range of the applied 
field, although their positions change continuously as a 
function of the field p]. Fig.|H|is a plot of the 2D ener- 
gyscape for a spherical particle with uniaxial anisotropy 
in the core, as before, but now with NSA on the surface. 
It is clear that the cubic-anisotropy features are also seen 
in the case of NSA, namely that i) the energy minima 
are not along the directions {9n = 0, tt) of the core easy 
axis; these directions having become local maxima, as 
discussed in the case of TSA, and ii) there is a clear de- 
pendence on the azimuthal angle ipn. In addition, more 
extensive calculations jllj have shown that similar fea- 
tures are also observed for other crystal structures (fee, 
bcc, etc.). 

These results agree and complement those of Ref. jj] , 
where it was shown that for the TSA and NSA there ex- 
ists a (different) critical value of the surface anisotropy 
constant that separates i) the OSP Stoner-Wohlfarth 
(SW) regime of coherent switching and ii) the MSP 
regime where the strong spin non-coUinearities invalidate 
the coherent mechanism, and the particle can no longer 
be modeled by an effective OSP. Obviously, for very small 
surface anisotropy the cubic contribution becomes negli- 
gible [see the first panel in Fig.01 kg = 0.1] and the OSP 
SW model provides a good approximation to the MSP. 
Accordingly, some experimental macroscopic estimations 
of the surface anisotropy constant yield, e.g., for cobalt 
kg ~ 0.1 foi' i^'Oii — 0.06 ^-iid for maghemite 
particles kg ~ 0.04 However, one should not forget 
that this effective constant depends on the particle's size, 
among other parameters such as the material composi- 
tion, and for, e.g. a diameter of 2 nm we may expect 
stronger anisotropies. 



IV. SURFACE EFFECTS ON THE DYNAMICS 

In the particular (and typical) situation, considered 
above, where the exchange interaction is much stronger 
than the anisotropy energy we find that the magneti- 
zation can be represented as the spatially homogeneous 
"global" magnetization plus a small inhomogeneous con- 
tribution that we calculated analytically and numerically 
0, 0, Q ■ The latter is induced by surface anisotropy and 
it is maximal near the surface but can extend deeply 
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FIG. 4: 2D energyscape for a spherical particle of TSA with varying constant ks- The symbols are results of the numerical 
calculations for the MSP, and the full lines are fits using Eq. HOII . (pn = and h = 0. 




FIG. 5; Plot of (-^ cos ip, ^ sin (p, E) for the same parameters 
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FIG. 6: Energy as a function of the polar angle On for a 
spherical particle with uniaxial anisotropy in the core and 
Neel anisotropy on the surface. (f>n = 0, 7r/4 and h = Q. 



into the body of the particle. It describes the adjust- 
ment of the magnetization to the conditions at the sur- 
face by minimizing the total energy with fixed direction 
of the global magnetization. As a result, we obtain the ef- 
fective particle's energy that depends on the orientation 
of its net magnetization and arises because of surface 
anisotropy. This contribution is of second order in the 
surface anisotropy and it adds to other terms, such as the 
bulk anisotropy and the first-order contribution from the 



surface anisotropy, which disappear in samples of cubic 
or spherical shape. These contributions to the energy of a 
magnetic nanoparticle are crucial to its dynamical behav- 
ior, in particular, in the ferromagnetic resonance (FMR). 
Accurately taking all of them into account should make 
it possible to determine the bulk and surface anisotropics 
from the experimental data. An interesting problem is 
the dynamical aspect of the magnetization adjustment 
mentioned above. If the anisotropics are much smaller 
than exchange interaction, the exchange-driven adjust- 
ment is much faster than the global precession of the 
magnetization induced by the anisotropy. Then, these 
adjustment modes behave adiabatically at low frequen- 
cies and the effective OSP energy is a good approxima- 
tion. For materials with a very strong surface anisotropy 
such separation of dynamical scales is no longer valid, 
and the dynamics of such nanoparticles becomes an es- 
sentially many-body process. FMR experiments on mag- 
netic nanoparticles should allow to estimate the values 
of the surface anisotropy and detect different regimes of 
their dynamical behavior. 

The effective OSP energy H10|l could be used in in- 
vestigating the thermally activated reversal of a small 
nanoparticle. Indeed, this effective energy, together with 
allows us to include surface effects, though in an 
approximate way, while offering a considerable simplifi- 
cation over the initial many-spin system. However, one 
should note that the effective potential energy contains a 
quartic term in the net magnetization components, which 
renders the analysis of the energyscape somewhat more 
involved but still tractable llSll. as opposed to the case of 
a many-spin particle. In Ref.ll^it was shown that for the 
case of mixed uniaxial and cubic anisotropics, there are 
two different relaxation rates, one for the parallel compo- 
nent of the magnetization and the other for the compo- 
nent perpendicular to the core uniaxial anisotropy axis. 
Using Eq. Hll|l one finds, for instance, that the parallel re- 
laxation rate increases with increasing surface anisotropy 
while the perpendicular component has the opposite be- 
haviour. This issue requires a more thorough investiga- 
tion. Nevertheless, for small enough surface anisotropy 
for which the effective energy H1U|I holds best, the minima 
are mainly defined by the uniaxial anisotropy. The effect 
of the cubic-anisotropy contribution is then to modify the 



6 



loci of the saddle points and their number, and thereby 
the calculation of the relaxation rates becomes more in- 
volved [13. 

V. CONCLUSION 

We have studied the properties of magnetic nanopar- 
ticles with the many-spin approach including the inter- 
nal structure and physical parameters such as core and 
surface anisotropy and exchange interactions. We have 
shown that upon varying the surface anisotropy con- 
stant there appear three regimes [see, e.g.. Fig. |2]: i) 
When this constant is very small with respect to the 
exchange coupling all observables can be scaled recov- 
ered within the Stoner-Wohlfarth macroscopic approach, 
ii) For small and intermediate values, the energy, and 
thereby the related quantities, of the many-spin parti- 
cle can be modeled by an effective expression contain- 
ing a surface-renormalized uniaxial contribution and a 
4*''-order contribution that is induced by the surface 
anisotropy. iii) For large values of the surface anisotropy 
constant the spin non-coUinearities are so strong that the 
magnetic state of a nanoparticle can no longer be repre- 
sented by a net magnetization. In this case, the switching 
of the particle's magnetization proceeds cluster-wise. 

In the intermediate regime, we have also discussed the 
possibility of studying the effect of surface anisotropy on 
the relaxation rate and speculated how to determine the 
nature and intensity of the surface anisotropy using the 
technique of ferromagnetic resonance. 

However, the present models on which all calculations 
are based suffers from some deficiencies: For instance, 
it is assumed that the crystal structure on the surface 
is the same as in the core with the same atomic lattice 
parameters. This cannot be wholly true considering the 



possibility of surface reconstructions. Of course, these 
models do include apices, edges and facets, and the pos- 
sibility of taking the exchange coupling on the surface 
as different from that in the core, or that between the 
core and the surface. However, there is no guarantee 
that it will become experimentally possible in the near 
future to estimate the atomic positions and lattice pa- 
rameters and may be theoretically possible to perform 
crystal field calculations, and eventually check these as- 
sumptions. Likewise, the magneto-crystalline anisotropy 
constant and exchange coupling in the core are taken as 
those in the bulk underlying material. However, we have 
shown that the core of a nanoparticle does not enjoy the 
properties of the underlying bulk material. For a given 
material such parameters could vary with the radius of 
the particle. In addition, the intensity and nature of sur- 
face anisotropy constitute a real challenge. There are 
many estimates but no firm understanding is achieved 
yet as to how surface anisotropy stems from the atomic 
structure in nanoparticles of reasonable size. 

For these interrogations to receive answers experimen- 
talists will have to devise ever more sensitive equipments 
and measurement techniques in order to probe the in- 
trinsic properties of nanomagnets and have direct access 
to the related physical observables. This, of course, will 
require new strategies for further isolating the nanomag- 
nets and rid them off the influence of their hosting matri- 
ces and mutual interactions. On the other hand, under- 
standing the influence on the particles static and dynamic 
behavior of the surrounding matrix and the inter-particle 
interactions, is of paramount importance to efhcient prac- 
tical applications. 

New experiments on adequately prepared particles are 
needed to check the validity of these results which, if con- 
firmed, should help us better understand the dynamics 
of nanoparticles beyond the macrospin approximation. 
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